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Abstract — The Sentinel-2 mission is a constellation of two satellites, developed by ESA on 
behalf of the joint ESA/European Commission initiative GMES (Global Monitoring for 
Environment and Security). In addition to monitoring land cover and disasters, it is useful 
in the assessment of the quality of terrestrial and marine waters with biophysical and 
chemical-physical parameters such as chlorophyll "a" and total suspended matter 
concentrations. The twin Sentinel-2A and Sentinel-2B satellites are equipped with a MSI 
(MultiSpectral Instrument) optical sensor that acquires in 13 spectral bands, with global 
coverage, a revisit time of 5 days and a spatial resolution of up to 10 meters on the ground. 

In this study, the estimates of chlorophyll "a" and of dispersion of sediments in the 
sea, derived from Sentinel-2 images were compared with in situ data acquired by means of 
a multiparametric probe in the monitoring stations that ARTA Abruzzo monthly checks 
along the Adriatic coast. The ultimate goal is to optimize parameters and algorithms to be 
able to derive concentration maps of chlorophyll and suspended solids from satellite, taking 
advantage of the high frequency of revisit time and high spatial resolution of the satellite 
acquisitions. This information is particularly relevant for aquaculture activities, for the 
monitoring of water quality and for the analyses of sedimentary processes. 

Sentinel-2 level 1C images were processed using the C2RCC (Case-2 Regional Coast 
Colour) processor in ESA SNAP software, setting the required parameters of the algorithm 
according to the in situ sampling measurements, in 4 acquisition dates during the year 2018. 

The comparisons, both at the sampling point and in the 100 meters radius buffer, 
showed statistically significant results: for chlorophyll “a”, linear regression models have 
r? = 0.82 in situ, r? = 0.81 in the 100 m buffer, respectively; r? = 0.88 in situ and r° = 0.91 
for TSM. 

This preliminary study has provided encouraging results with only four sampling 
dates. In future ongoing work, the increase in the number of sampling dates should further 
improve the accuracy of the estimates; we aim also to compare other satellite products (e.g. 
Sentinel-3 OLCI) and test other algorithms (ACOLITE, BOMBER, etc.). All these aspects 
should lead to obtain a most accurate regional algorithm and an efficient process, for the 
production of maps and indicators useful in the monitoring of water quality. 
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Introduzione 


Il monitoraggio delle acque costiere è richiesto da direttive europee [7, 8, 9, 23] e 
da regolamenti nazionali sulle attività produttive nelle aree litoranee: impianti di 
trattamento delle acque, attività di estrazione (ad es. petrolio e gas), acquacoltura, etc. Il 
raggiungimento del “buono” stato di qualità dell'ambiente marino è quindi un obiettivo 
comunitario e nazionale, parte della strategia di salute globale “One Health”, condizionata 
dalle attività antropiche e da cambiamenti climatici. 

La necessità di dati sulla qualità delle acque (in senso lato) è stata finora 
soddisfatta principalmente tramite campionamenti e misurazioni di parametri chimico-fisici 
effettuati in mare [22]. Questo rende le campagne di acquisizione dati onerose in termini di 
risorse impiegate e permette campionamenti limitati nel tempo e nello spazio, mentre 
normalmente le aree marine da monitorare sono ampie ed è necessario un monitoraggio 
costante e frequente. Campagne di campionamenti e sistemi di monitoraggio per la qualità 
dell’acqua e la salubrità dei prodotti ittici in aree costiere sono utilizzati in acquacoltura 
[11, 24, 25], nel trattamento delle acque reflue, nelle attività di dragaggio, nella protezione 
delle Aree Marine protette, in tutte le attività commerciali e turistiche che insistono e 
convivono sulle zone costiere [10, 21]. I parametri rilevati e misurati, comuni a queste 
finalità, sono, principalmente, clorofilla “a”, trasparenza delle acque, torbidità, eventuale 
presenza di bloom algali e solidi sospesi. 

Inoltre, l’Osservazione della Terra è una delle principali conquiste scientifico- 
tecnologiche degli ultimi decenni, poiché fornisce dense serie temporali di dati (con 
passaggi frequenti dei satelliti) su vaste aree geografiche (big swath, ossia grande superficie 
acquisita a terra) rilevando la radiazione riflessa dalle superfici in ampi intervalli dello 
spettro elettromagnetico. 

Molti sono i prodotti di remote sensing dedicati al monitoraggio delle acque marine 
(es. Copernicus Marine Environment Monitoring Service - CMEMS, MODIS, Sentinel-3 
OLCI). La loro media risoluzione spaziale (dimensione a terra della porzione di territorio 
rappresentata da un pixel) varia generalmente da 300 metri a 1 km, o dimensioni maggiori. 
Sebbene i sensori multispettrali a media risoluzione spaziale montati sulle piattaforme 
satellitari abbiano una configurazione radiometrica più adatta alla stima di parametri 
descrittivi della colonna d’acqua (es. Chl-a e TSM), la loro risoluzione spaziale rende 
difficilmente utilizzabile i loro prodotti in applicazioni di monitoraggio in prossimità delle 
coste, per le quali si rende necessaria l’identificazione di pattern spaziali anche a piccola scala. 

I satelliti Sentinel, realizzati nell’ambito del programma comunitario Copernicus [20], 
contribuiscono al sistema globale di monitoraggio per l’osservazione della Terra (Global 
Earth Observation System of System - GEOSS); il programma è coordinato dalla 
Commissione Europea con la compartecipazione dell’ Agenzia Spaziale Europea (ESA). 
L’operatività è assicurata da un insieme di sei missioni (da Sentinel-1 a Sentinel-6) 
ciascuna dedicata a specifiche applicazioni. I satelliti Sentinel-2, nello specifico, oltre al 
monitoraggio della copertura del suolo e dei disastri, acquisiscono in prossimità della costa 
e possono essere impiegate per il monitoraggio delle acque terrestri e marine e per la 
mappatura delle variabili chimico-fisiche [16, 26], come la clorofilla “a” [24] e le 
concentrazioni di solidi sospesi totali (TSM) [2, 4, 5]. 

In questo studio, i parametri caratterizzanti la colonna d’acqua, clorofilla “a” e 
solidi sospesi totali, sono stati stimati da immagini Sentinel-2 MSI ed è stata confrontata 
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con dati in situ acquisiti in mare da sonda multiparametrica, lungo la costa abruzzese nel 
Mare Adriatico. L’obiettivo finale è quello di ottimizzare parametri ed algoritmi per poter 
derivare mappe di concentrazione di clorofilla “a” e solidi sospesi da satellite, sfruttando la 
loro elevata frequenza temporale e l’alta risoluzione spaziale. La generazione di tali 
informazioni potrebbe essere di rilevante supporto all’attività di acquacoltura, per il 


monitoraggio della qualità delle acque e per l’analisi dei processi sedimentari. 


Materiali e Metodi 


Area di studio 


L’area di studio è lo specchio di Mare Adriatico prospiciente le coste della Regione 
Abruzzo. In questo tratto, lungo circa 120 km, si versano numerosi corsi d’acqua nelle cui 
pianure alluvionali di fondovalle sono presenti molte delle principali zone industriali ed 
agricole della regione. Le dinamiche di dispersione nel Mar Adriatico sono influenzate 
dalle correnti, dalle acque provenienti dal Mar Ionio, dalla portata del fiume Po e 
soprattutto dalla piovosità che provoca il maggior apporto di inquinanti in mare [6]. 


Dati rilevati in mare 


Il piano di monitoraggio istituzionale che l ARTA Abruzzo conduce sul controllo 
delle acque marino-costiere, è espletato sulla rete regionale, come indicato dalla convenzione 
con il Settore Opere Marittime della Regione Abruzzo, su 14 stazioni di campionamento 
distribuite su 7 transetti perpendicolari alla costa, ai fini della determinazione dello stato di 
qualità ambientale dei corpi idrici superficiali come previsto dal Direttiva 2000/60/CE. In 
ciascuna stazione sono state effettuate rilevazioni chimico-fisiche sulla colonna d’acqua: 
temperatura, salinità, conduttività, clorofilla “a”, torbidità, ossigeno disciolto e pH. 

Questi dati sono stati rilevati tramite sonda multiparametrica IDRONAUT Ocean 
Seven 316 Plus equipaggiata con fluorimetro e torbidimetro Seapoint Sensors, inc. 


Dati satellitari 


Le concentrazioni di clorofilla (CHL) e di solidi sospesi totali (TSM) in superficie, 
misurati in situ, sono state appaiate a stime derivate dalle misure radiometriche degli 
strumenti a bordo dei satelliti. 

Il dataset di immagini satellitari utilizzati è quello della missione Copernicus 
Sentinel-2, con i due satelliti gemelli, Sentinel-2A e Sentinel-2B, che assicurano dati con 
tempi di rivisitazione (sullo stesso punto) e copertura globale di 5 giorni (ciascun satellite 
ripassa sullo stesso punto ogni 10 giorni) [16]. Questi satelliti sono dotati di un sensore 
ottico MSI (MultiSpectral Instrument) che acquisisce in 13 bande spettrali, tra 1443 nm e i 
2190 nm, con una risoluzione spaziale tra 10 e 60 metri al suolo. 

Per ciascun prelievo in mare è stato acquisito il passaggio dei satelliti Sentinel-2A o 
Sentinel-2B avvenuti lo stesso giorno (con una tolleranza di qualche ora tra il passaggio del 
satellite e la misurazione in mare) e senza copertura nuvolosa. 


Processamento dei dati satellitari 


I dati Sentinel-2 MSI sono stati elaborati in tre fasi: acquisizione, elaborazione e 
analisi statistica. 
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L’acquisizione è avvenuta tramite un batch script Python che sfrutta, previa 
registrazione, il servizio Copernicus Open Access Hub [18], grazie al pacchetto ufficiale 
Sentinelsat, che scarica il dato Sentinel-2 nel formato compresso SAFE (Standard Archive 
Format for Europe). 

I dati presi in esame (livello 1C) ricadono nelle granule T33TVG e T33TVH della 
griglia di distribuzione, appartenenti alle orbite 79 e 122 dei satelliti. 

Il livello 1C fornisce immagini ortorettificate di dimensione 100 km x 100 km, nella 
proiezione Universal Transverse Mercator/World Geodetic System 1984 (UTM/WGS84) e 
contiene la riflettanza Top of Atmosphere, ossia non corretto dall’effetto dell’atmosfera. 

La fase di elaborazione, effettuata con il software SNAP v 6.0, distribuito da ESA, 
ha previsto il ricampionamento (o resampling) di tutte le bande alla risoluzione spaziale di 
10 m, e l’estrazione delle bande dedicate di clorofilla “a” e solidi sospesi, prodotte dal 
processore di colori oceanici multi-missione C2RCC S2-MSI (Case-2 Regional 
CoastColor) [3], specifico per le Sentinel-2 e risultato tra i migliori stimatori di questi 
parametri in studi simili [1, 17]. È stato utilizzato l'algoritmo C2RCC-Nets, parametrizzato 
con i valori di salinità e temperatura ottenuti dai campionamenti in situ (tabella 1), e dalla 
pressione atmosferica relativa all’ Aeroporto di Pescara [19] ricalcolata all’altezza del 
livello del mare. Dalle bande ottenute di Chlorophyll concentration (conc_chl) e total 
suspended matter dry weight concentration (conc_tsm), sono stati estratti sia i valori nelle 
stazioni di campionamento, sia la media dei valori in un buffer di 100 metri di raggio 
intorno a ciascuna stazione, per ciascuna data di acquisizione. Per questa elaborazione è 
stato utilizzato il software ESRI® ArcMap 10.5. 

Sui dati estratti sono state effettuate analisi di regressione lineare per verificare 
l’esistenza di una relazione tra il valore dei parametri CHL e TSM stimati da rilevazioni 
satellitari e i valori degli stessi parametri campionati in mare. L’analisi statistica è stata 
effettuata mediante l’utilizzo del software R Core Team (2019). 


Risultati 


La figura 1 mostra la localizzazione delle 14 stazioni di campionamento dislocate 
lungo la costa abruzzese del Mar Adriatico, e sovrapposte, le delimitazioni delle immagini 
Sentinel-2A e Sentinel-2B (granule T33TVG e T33TVH). 

Considerate le date di campionamento ed i passaggi dei due satelliti, sono state 
processate 4 immagini satellitari in 4 diverse giornate di campionamento in mare, per un 
totale di 18 osservazioni in 12 punti (Tabella 1). Tre osservazioni sono state escluse dalle 
regressioni poiché outlier nel dataset. Inoltre, per tener conto dei potenziali minimi 
scostamenti delle coordinate dei punti di campionamento, della sovrapposizione delle 
granule e della variabilità locale dei parametri misurati, è stato considerato sia il valore ne 
punto sia la media dei pixel in un intorno di 100 metri di raggio attorno a ciascun punto. 

Ai dati estratti dall’elaborazione delle immagini satellitari, sia nel punto di 
campionamento, sia nel suo intorno di 100 metri di raggio, è stato applicato un modello di 
regressione lineare, che ha evidenziato risultati statisticamente significativi sia per la CHL 
(r? = 0,82 nel punto, figura 2; r? = 0,81 nel buffer di 100 metri, figura 3) sia per i TSM 
(r? = 0,88 nel punto figura 4, r?° = 0,91, figura 5). 
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Applicando tali regressioni (media nel buffer a 100 metri) alle quattro immagini 
considerate, sono state ottenute le mappe di CHL (Figura 6) e TSM (Figura 7). 


Punti di campionamento 
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Figura | - Localizzazione dell’area di studio e delle stazioni di campionamento. I riquadri 
mostrano l’extent delle immagini satellitari acquisite dai sensori a bordo dei satelliti 
Sentinel-2A e Sentinel-2B. 

Figure 1 - Study area (coast facing Abruzzo region in the Adriatic Sea) and location of the 
sampling points. The squares are the extent of satellite imagery of Sentinel-2A e Sentinel-2B. 


Tabella 1 - Campionamento in situ: date, stazioni di campionamento e misure rilevate, 
utilizzate come parametri del modello. 


Table 1 - In situ sampling: dates, sampling points and in situ collected values used in the 
models. 


Data di Stazioni di Salinità | Temperatura Pressione Satellite 
campionamento | campionamento |media (PSU)| media (°C) |atmosferica (hPa) |corrispondente 
ORO7, OR09, VA10, - 
03-01-2018 VA12, SS01, SS02 36,5 9,2 1016,5 Sentinel-2A 
09-03-2018 PI16, PI18 32,3 9,3 1019,5 Sentinel-2B 
ORO7, OR09, VA10, ; 
20-06-2018 VA12. SS01, S02 35,2 24,8 1020,5 Sentinel-2B 
12-07-2018 pane STILL 34,8 21,1 1016,5 Sentinel-2A 
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Figura 2 - Modello di regressione lineare tra dati campionati in mare (in 
situ) e le concentrazioni di clorofilla stimate a partire dai dati satellitari 
(Satellite) sul punto di campionamento (mg/m). 

Figure 2 - Linear regression model between sampling collected in the sea 
(in situ) and the estimated concentration in the same point of chlorophyll 
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Figura 3 - Modello di regressione lineare tra dati campionati in mare (in situ) e 
le concentrazioni di clorofilla stimate a partire dai dati satellitari (Satellite 
100 m) nell’intorno di 100 m dal punto di campionamento (mg/m?). 
Figure 3 - Linear regression model between sampling collected in the sea (in 
situ) and the average estimated concentration in a 100 m - buffer around the 
point of chlorophyll “a” derived from satellite imagery (Satellite) (mg/m’). 
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Figura 4 - Modello di regressione lineare tra dati campionati in mare (in situ) con 
fluorimetro (in FTU) e le concentrazioni di solidi sospesi stimate con l’algoritmo 
C2RCC dai dati satellitari (Satellite) sul punto di campionamento (in g/m’). 
Figure 4 - Linear regression model between sampling collected in the sea (in 
situ) measured with a fluorimeter (in FTU) and the concentration in the same 
point of Total Suspended Matter (TSM) estimated with the C2RCC algorithm 
from satellite imagery (Satellite) (in g/m’). 
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Figura 5 - Modello di regressione lineare tra dati campionati in mare (in situ) 
con fluorimetro (in FTU) e le concentrazioni di solidi sospesi stimate con 
l’algoritmo C2RCC dai dati satellitari (Satellite 100 m) nell’intorno di 100 m 
dal punto di campionamento (in g/m?). 

Figure 5 - Linear regression model between sampling collected in the sea (in 
situ) measured with a fluorimeter (in FTU) and the average concentration in 
a 100 m-buffer around the point of Total Suspended Matter (TSM) estimated 
with the C2RCC algorithm from satellite imagery (Satellite 100 m) (in g/m’). 
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Figura 6 - Mappe di concentrazione della clorofilla (Chl-a) della costa abruzzese elaborate dalle 
immagini Sentinel-2 con il processore C2RCC per le date 03/01, 09/03, 20/06 e 12/07 del 2018, 
corrispondenti alle granule T33TVG (nelle mappe di sinistra) e T33TVH (nelle mappe a destra). 
Le mappe sono ottenute applicando la regressione lineare al buffer a 100 metri (figura 3). 
Figure 6 - Concentration maps of chlorophyll (Chl-a) along Abruzzo coast in the Adriatic 
sea, as elaborated with the C2RCC processor from Sentinel-2 imagery of the dates: 
03/01/2018, 09/03/2018, 20/06/2018 and 12/07/2018. The values in the maps are obtained 
by applying the linear regression in the 100 m — buffer (figure 3). The images on the left 
side are the T33TVG granule, the images on the right are the T33TVH granule. 


Discussione 


Questo studio preliminare, con sole 4 campagne di acquisizione nel 2018, 
simultanee ad altrettanti passaggi dei satelliti Sentinel-2 (A o B), ha fornito risultati 
incoraggianti per la stima da dati telerilevati di alcuni parametri bio-geofisici dell’acqua 
lungo la costa abruzzese. Il settaggio dei parametri di stima dell’algoritmo deve ancora 
essere migliorato, soprattutto per quanto riguarda le performances in termini di accuratezza. 
Al momento, i pochi dati disponibili mostrano comunque l’esistenza di una correlazione, 
incoraggiandoci nella ricerca della migliore calibrazione possibile. Lo studio, tutt'ora in corso, 
prevede l’incremento della numerosità dei dati raccolti in situ per rendere le stime più 
accurate, la comparazione con altri prodotti satellitari (es. Sentinel-3 OLCI [12] e LANDSAT8 
OLI [15]) e l’utilizzo di altri algoritmi (ACOLITE, BOMBER, etc.) per valutare il processo più 
accurato, efficiente (tempi di elaborazione, complessità e regionalizzazione tramite procedure 
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Figura 7 - Mappe di concentrazione dei sedimenti di solidi sospesi (TSM) della costa 
abruzzese elaborate dalle immagini Sentinel-2 con il processore C2RCC per le date 03/01, 
09/03, 20/06 e 12/07 del 2018 con i punti di campionamento in situ, corrispondenti alle 
granule T33TVG e T33TVH. 

Figure 7 - Concentration maps of Total Suspended Matter (TSM) along Abruzzo coast in 
the Adriatic sea, as elaborated with the C2RCC processor from Sentinel-2 imagery of the 
dates: 03/01/2018, 09/03/2018, 20/06/2018 and 12/07/2018. The values in the maps are 
obtained by applying the linear regression in the 100 m — buffer (figure 5). The images on 
the left side are the T33TVG granule, the images on the right are the T33TVH granule. 


di calibrazione), in grado di correggere effetti relativi alle condizioni presenti al momento 
dell’acquisizione delle immagini (es. swn-glint) e capace di produrre mappe ed indicatori per 
il monitoraggio della qualità delle acque. 

Le unità di misura non sono sempre corrispondenti, dato che generalmente le 
misure in situ per la caratterizzazione della qualità delle acque vengono fatte con strumenti 
che misurano fluorescenza e torbidità, mentre gli algoritmi da satellite stimano Chl e TSM. 
Campagne ad hoc saranno necessarie per misurare in situ il dato da torbidimetro, e 
contestualmente effettuare prelievi di acque superficiali da analizzare poi in laboratorio. 
Così facendo avremo la possibilità di confrontare il dato satellitare con quello reale 
utilizzando la stessa unità di misura prevista dai metodi analitici per le acque. 

L’obiettivo è di regionalizzare gli algoritmi di stima per le acque della costa 
abruzzese, con risultati coerenti e affidabili, che potranno in futuro essere utilizzati per 
ottimizzare gli strumenti di monitoraggio riducendo anche l’impiego delle risorse in mare. 
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A regime si potranno avere mappe multitemporali di clorofilla “a” e solidi sospesi totali su 
vaste porzioni di superfici [22], finalizzati al monitoraggio dei loro pattern spaziali di 
distribuzione [14, 13]. La rilevanza nell’utilizzo dei dati Sentinel-2 é data dalla risoluzione 
spaziale (10 metri), compatibile con la scala di osservazione necessaria alla valutazione e al 
monitoraggio degli effetti nel mare delle attivita antropiche. 
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